Ultrastable Optical Clock with Neutral Atoms in an Engineered Light Shift Trap 
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An ultrastable optical clock based on neutral atoms trapped in an optical lattice is proposed. 
Complete control over the light shift is achieved by employing the 5s 2 ^So — > 5s5p 3 Po transition of 
87 Sr atoms as a "clock transition". Calculations of ac multipole polarizabilities and dipole hyper- 
polarizabilities for the clock transition indicate that the contribution of the higher-order light shifts 
can be reduced to less than 1 mHz, allowing for a projected accuracy of better than 10~ 17 . 

PACS numbers: 32.80.Pj, 42.50.Vk, 31.10.+Z, 31.15.Ar, 31.15.Md, 32.70.Cs 



O 
-I— > 

ctf ; 

O ■ 

43 ' 

Oh: 



> 
m 

o 

On 
O 
CO 

o 

o 

43 

Or 



X 



Careful elimination of perturbations on electronic 
states and of motional effects has been considered as a 
prerequisite for realizing an atom frequency standard . 
A single ion trapped in an RF quadrupole field is one 
of the ideal systems that satisfy these requirements 0, 
as the trap prepares a quantum absorber completely at 
rest in free space for an extended time and its electric 
field vanishes at the center of the trap. Employing this 
scheme, quantum projection noise (QPN) limited spec- 
troscopy [3j has been performed with an expected accu- 
racy of 10" 18 [10. 

Despite its anticipated high accuracy, the stability of 
the single-ion based optical clock is severely limited by 
QPN; long averaging times are required to meet its ulti- 
mate accuracy |5j . The measure of the fractional instabil- 
ity is provided by the Allan variance, <J y {r) = ^—j===. 

Assuming the transition line Q ~ 1.6 x 10 14 Q and a cy- 
cle time of T m f» 0.1s, 4 x 10 7 measurement cycles are 
required for a single quantum absorber (N — 1) to reach 
<J y {r) = 10 -18 , corresponding to a total averaging time r 
of a few months. For further increase of the stability, the 
averaging time increases quadratically and will become 
inordinately long. 

One may think of increasing the number of quantum 
absorbers N as employed in neutral atom based optical 
standards In this case, however, the atom-laser 

interaction time sets an upper bound for the Q-factor 
since an atom cloud in free space expands with finite ve- 
locity and is strongly accelerated by the gravity during 
the measurement. Hence the highest line Q « 10 12 ob- 
tained for neutral atoms is 2 orders of magnitude smaller 
than that of a trapped ion. Furthermore, it has been 
pointed out that residual Doppler shifts arising from an 
imperfect wavefront of the probe beam and atom-atom 
collisions during the measurement affect its ultimate ac- 
curacy 



In this Letter, we discuss the feasibility of an "opti- 
cal lattice clock" which utilizes millions of neutral 
atoms separately confined in an optical lattice [10] that 
is designed to adjust the dipole polarizabilities cxei for 
the probed electronic states in order to cancel light field 
perturbations on the measured spectrum |ll| . In striking 
contrast with conventional approaches toward frequency 
standards [l|, the proposed scheme interrogates atoms 
while they are strongly perturbed by an external field. 
We will show that this perturbation can be canceled out 
to below 10 -17 by carefully designing the light shift po- 
tentials. This scheme permits an exceptionally low in- 
stability of Cj,(t) « 10~ 18 with an interrogation time of 
only r = 1 s, which may open up new applications of 
ultra precise metrology, such as the search for the time 
variation of fundamental constants |l2| and the real time 
monitoring of the gravitational frequency shift. 

Figure ^ illustrates the proposed scheme. Sub- 
wavelength confinement provided by the optical lattice 
localizes atoms in the Lamb-Dicke regime (LDR) [l3j . 
where the first order Doppler shift as well as the pho- 
ton recoil shift disappears 0, 0] and the second order 
Doppler shift can be made negligibly small by sideband- 
cooling atoms down to the vibrational ground state 
0, fl~5j | . In addition, a 3-dimensional lattice with less 
than unity occu pat ion could reduce the collisional fre- 
quency shifts [a, lla] • Therefore this scheme simulates a 
system where millions of single-ion clocks operate simul- 
taneously. 

The transition frequency v of atoms exposed to the 
lattice electric field of £ is described as, 

hv = h»V> - iAa(e, co)£ 2 - ^A 7 (e, c)£ 4 - . . . , (1) 

where is the transition frequency between the un- 
perturbed atomic states, Aa(e,w) and A7(e,u;) are the 
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FIG. 1: Simplified optical coupling scheme for 87 Sr. (a) In 
the limit of large detunings Si of the coupling laser compared 
to the hyperfine splittings <5hfs, the squared transition dipole 
moment of the upper J manifold can be simply added up, re- 
sulting in a quasi-scalar light-shift, (b) 3D optical lattice pro- 
vides Lamb-Dicke confinement while it prevents atom-atom 
interactions. 

difference between the ac polarizabilities and hyperpo- 
larizabilities of the upper and lower states, which in the 
general case depends both on the light wave frequency ui 
and on the unit polarization vector e. Higher order cor- 
rections are included in the hyperpolarizability 7(e, ui) 
and in the higher-order multipole corrections to the po- 
larizability (magnetic dipole Ml and electric quadrupole 
E2 terms in addition to the electric dipole polarizability 

a(e,ui) = a E i{e,u)) + a M i{e,u>) + a B2 (e,cj). (2) 

By canceling out the polarizabilities of the upper and 
lower states to set Aa(e,u>) — 0, the observed atomic 
transition frequency will be equal to the unperturbed 
transition frequency independent of the laser intensity 
I oc \£\ 2 , as long as higher order corrections 0(£ 4 ) are 
negligible. 

We first discuss the cancellation of the electric dipole 
polarizabilities aEi(e,w) that dominate the light shift 
perturbations. In order to provide complete control over 
perturbations, parameters are preferably specified by the 
frequency u> that is the most accurately measurable value 
in experimental physics. Other parameters such as the 
light polarization e should have less influence on the sys- 
tem: Our strategy is to employ the J = state that 
exhibits a scalar light shift Q. Here we propose the 
5s 2 -"So — * 5s5p 3 Po forbidden transition (Ao = 698 nm) 
of 87 Sr with nuclear spin I = 9/2 as the "clock" tran- 
sition 9], in which we take advantage of the hyperfine 
mixing of the 3 P (F = 9/2) state with the ^Pi states 
[l7l | to gain a finite lifetime of Tq 1 = 150 s. This tran- 
sition rate turned out to be ideal to perform ultra nar- 
row spectroscopy, as the blackbody radiation at 300 K 
quenches the 3 P metastable states of strontium and gives 



an effective lifetime of ~ 100 s [lSg. Experimentally, how- 
ever, the attainable line Q-factor will be limited by the 
collision- limited lifetime of 10 s at the vacuum pressure 
of 10~ 10 torr and/or by the mean photon-scattering time 
(~ 10 s) of atoms in a far-off-resonant light trap. 

Figure |21 shows the light shift for the 1 5o and 3 Po states 
as a function of the trapping laser wavelength with an 
intensity of I = 10kW/cm 2 . The calculation is per- 
formed by summing up the light-shift contributions with 
electronic states up to n = 11 orbits [Jjj> m which wc 
employed new values of the dipole moments determined 
in recent experiments [lij to find the intersection wave- 
length at Al ~ 800 nm. The light shifts around the inter- 
section are mainly determined by the states indicated in 
Fig.^a): The light shift of the 3 Po state can be tuned ar- 
bitrarily in the near infrared range, being "pushed down- 
ward" by 3 5i and "pulled upward" by the 3 P>i state, while 
that of the x So state monotonically decreases toward the 
dc Stark shift. This tuning mechanism, therefore, can be 
similarly applied to heavier alkaline-earth (like) atoms. 
At Al the light shift v ac changes with the trapping laser 
frequency u> as d^ ac /dcj = —3.6 x 10~ 10 and —1.3 x 10~ 9 
for the and 3 Po state, respectively. This precision en- 
hancement of more than a factor of 10 9 allows to control 
the light shift well within 1 mHz by defining the cou- 
pling laser within 1 MHz or 10~ 8 precision for the optical 
frequency, which can be easily accomplished by conven- 
tional optical frequency synthesis techniques. Using the 
same set of dipole moments as used in the above cal- 
culation, the blackbody frequency shift in the clock 
transition is estimated to be —2.4 Hz at T = 293 K: The 
uncertainty can be reduced to the 10 mHz level by con- 
trolling the surrounding temperature variation AT < 0.5 
K. Alternatively, by operating the clock at T = 77 K, 
the blackbody shift is dramatically decreased to 10 mHz 
because of its T 4 dependence. 

Owing to the hyperfine interaction (HFI), the clock 
transition no longer consists of simple scalar states; there- 
fore the tensor light shift arises. The fraction of the ten- 
sor shift in the total light shift is roughly proportional 
to 5hfs/S, i.e., the ratio of the hyperfine splitting of the 
coupled electronic state <5hfs to the trapping laser detun- 
ing 5 to that state. For the 3 5i, 3 D%, and Pi state, 
as shown in Fig. da), the ratio Sms/S is 8.1 x 10~ 5 , 
5.4 x 10~ 6 , and 2.1 x 10 -7 , respectively. Therefore the 
tensor shift of the 3 Po state can be 2 orders of magnitude 
larger than that of the ^So state. In order to give an 
approximate estimate for this tensor shift, we calculated 
the polarization-dependent light shift in the 3 Po state, 
taking into account the electric-dipole coupling to the 
3 5i and 3 P>i hyperfine manifolds that are energy shifted 
due to HFI. Both of these manifolds provide half (i=s 65 
kHz) of the total light shift but dominate its tensor com- 
ponent: Although higher lying electronic states support 
another 50 % of the light shift, their contribution to the 
tensor shift can be less than 50 %, as the ratio S^/S 
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FIG. 2: Light shift as a function of the trapping laser wave- 
length for a laser intensity of I = 10kW/cm 2 . The solid and 
dashed lines show the light shifts for the 1 So and 3 Prj states, 
respectively, which intersect at Al ~ 800 nm. The inset illus- 
trates the insignificance of the polarization-dependent light 
shifts of the 3 Pn (F = 9/2) state by taking into account the 
dipole coupling to the 3 Si and 3 Di state at \l in presence of 
a magnetic field Bq — 30 mG. 



monotonically decreases for larger 8. The inset of Fig. |3 
shows the result with the light polarization parameter- 
ized as e = cos 9 e + i sin 9 e + (e± represent the unit 
vector for polarization), where a bias magnetic field of 
B = 30 mG is applied to reduce the Raman coherences 
among the Zeeman substates. The m = ±1/2 state can 
be best used for the "clock" transition, as it exhibits the 
smallest polarization-dependence of less than 1 Hz. By 
employing these states, one could control the light shift 
within 1 mHz by defining 9 within 1 mrad, even if one 
applied the linearly polarized trapping laser (9 = 7r/4) 
where the gradient dv ac /d9 (= 0.83 Hz/rad) is the largest. 

In order to provide an estimate for the higher-order 
field contributions 0{£ 4 ) described by ac h ype rpolariz- 
abilities, we have used the general theory [iol Ell E2I 
|H for calculating the light shift for the 5s 2 l S -> 
5s5p 3 Po transition on the basis of the Green's func- 
tion method in the Fues' model potential approximation 
|2~l| . These model potential calculations have reproduced 
Fig. |21 to within 3 % accuracy both for the amount of 
the ac Stark shifts as well as the intersection frequency 
lul(= 27rc/Ai), confirming the validity of these two inde- 
pendent approaches. 

To calculate the Ml and E2 contributions to the 
polarizability Eq. (J2J), the magnetic dipole and electric 
quadrupole atom-field interactions should be taken into 
account together with the electric dipole term in the am- 
plitude, 



V(r) = Vei+Vmi + V E2 , 



(3) 



of the interaction Hamiltonian H(r,t) — V(r)e~ luJt + 
V^(r)e luJt 21]. The magnetic dipole polarizability for 
the ground state equals zero, while for the excited state 



5s5p 3 Po its value is proportional to the squared fine- 
structure constant (a = 1/137), to the splitting of the 
triplet states E\q = E 5p 3p 1 — E 5p 3p and to the square of 
the wave functions overlapping integral. The quadrupole 
polarizabilities of both levels are on the order of (auf) 2 
and may also become considerable only in the closest 
vicinity of the resonance on the quadrupole-allowed tran- 
sition. Numerical estimates with the frequency col where 
Aoiei{ul) — gives ajvn ~ &E2 ~ 10~ 7 x o>ei for both 
levels. 

The hyperpolarizability j(e,uj) is calculated starting 
from a formal expression for the fourth-order quasi- 
energy in terms of the field-free wave functions 10) and 
the reduced quasi-energy Green's functions G |2ll l22|: 



c4 

AE^ = ~— 7 (e,w) = -({0\HGHGHGH\0}} 
+ ((0\HGH\0))((0\HG 2 H\0)), 



(4) 



where the double brackets indicate the integration over 
the spatial variables and averaging over time. In Eq. Q 
only the dipole term of Eq. is taken into account in 
the interaction Hamiltonian V'(r). 

After the time averaging and angular integration, 
the dipole dynamic polarizability and hyperpolarizability 
tensors may in general be resolved into 3 and 5 irreducible 
parts, a p (p = 0,1,2) and j q (q — 0,1,2,3,4), respec- 
tively, of which only scalar terms an and T n c ontribute 



in a state with the total momentum J = 



Cncoi 



and 



7 g are determined by linear combinations of frequency- 
dependent radial matrix elements of the second, third 
and fourth orders. 

It is to note first that the scalar parts of 70 (e, a;) are 
different for the linear and circular type of polarization, 
7o(w) 7^ 7o( w )j even for atoms in ^-states, while all the 
terms of the polarizability in Eq. J5| are independent of 
e for a state with J = 0. Secondly, the number and 
type of singularities for the hyperpolarizability exceeds 
that for the polarizability, and the contribution of these 
singularities also depend on the polarization of the light 
field. E.g., for linearly polarized radiation, there are two- 
photon singularities of 7q(^) on the J = states, while 
for the circular polarization such singularities cannot ap- 
pear. So, the hyperpolarizability of the ground state with 
two equivalent electrons 5s 2 in Sr may be written as 



:cri2i(w,2w,w) 



7» 
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SlOl + 7^-121 + 70"i2l(w, 2U), to) 

5 5 



2a (u>)S- 3 (uj) 



(5) 



where an(w) is the polarizability and S*_3 is the so-called 
frequency-dependent oscillator strengths moment. The 
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following notations for the radial matrix elements and 
their combinations were used above: 



cr ;i ; 2 i 3 (w, 0, w) + a hh i 3 ((j, 0, -w); 
^2/3(^1,^2,^3) 

(0|r 5i > 5 £ 2 r<£ a r|0>. (6) 



The arguments in Eqs. (6) indicate the frequency de- 
pendence in energy denominators [2l|. Here, gf is the 
radial Green's function in the subspace of the jumping 
electron's states with angular momentum I. 

In the vicinity of a one-photon resonance with fre- 
quency detuning \S\ <C to, the product of the second- 
order and third-order matrix elements in Eq. which 
corresponds to the term aoS-3 in Eq. JSJ, dominates 
the third-order poles of order 5~ 3 . The fourth-order ma- 
trix element in Eq. @ has the second-order poles (S~ 2 ). 
The two-photon (2uj) resonance singularity on the J = 
states in the radial matrix elements i?ioi(w, 2uj,uj) ap- 
pears only for 7 ( (w) and that on the J — 2 states in 
Ri2i(u!,2uj,uj) appears both for 7*(w) and 7 c (w). 

We used the analytic Sturm-series representation of 
the Green's function, corresponding to calculation of the 
infinite sums over the total atomic spectrum including 
the continuum. Finally, the radial integrals in Eq. © 
were presented in the form of absolutely converging se- 
ries, well suited for the numerical computations. The 
numerical results for 7^ (u) and ji p (w) at the intersec- 
tion frequency of ojl are 6.3 • 10 6 a.u. and 2.7 • 10 8 a.u., 
respectively, which give fourth-order ac Stark shifts of 



AE^/h: 



-5.3 • 1(T 5 Hz and AE^Jh « -2.3- lCT a Hz 
for the trapping laser intensity of 10kW/cm 2 . Therefore 
the contribution of the higher order light shifts is as small 
as 5 x 10~ 18 . Further elimination of this systematic er- 
ror can be achieved by extrapolating the trapping laser 
intensity to zero |l5l | in a quadratic way. 

In summary, we have discussed the feasibility of pre- 
cision spectroscopy of neutral atom ensembles confined 
in an optical lattice, by applying a light-shift cancella- 
tion technique on the 1 S , (F = 9/2, rap = ±1/2) — > 
3 F (F = 9/2, m F = ±1/2) clock transition of 87 Sr that 
has a negligibly small tensor shift and a suitable transi- 
tion moment to perform spectroscopy. Our calculation 
including the higher-order Stark shift confirmed that the 
measurement of the unperturbed atomic transition at 1 
mHz level is feasible, allowing for the projected accuracy 
of better than 10 -17 . Since this scheme is equivalent to 
millions of single-ion-clocks operated in parallel, excellent 
improvement of the stability ~ VlO 6 can be expected. 
The theory and method for the numerical calculation of 



the hyperpolarizability tensor for arbitrary alkaline-earth 
atoms will be discussed in a forthcoming paper. 
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Note Added: After submission of the paper, we have 
experimentally demonstrated the optical lattice clock and 
determined the intersection wavelength [2f|. 
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